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PARAMETRIC ANALYSIS OF ELECTROSTATIC 
DISPERSION OF A LIQUID 


Allan Sherman 
Auxiliary Propulsion Branch 


ABSTRACT 

This document presents a theoretical model for the electro- 
static dispersion of a fluid. The analysis results in equations for 
charge-to-mass ratio, velocity, and droplet radius as a function of 
flow-voltage parameters. Data from colloid rocket tests agree 
favorably with theory. The results are useful for colloid-rocket 
data correlation. 
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NOTATIONS 


A flow area of capillary (m 2 ) 

C wave velocity 
D diameter of capillary (m) 

f(e, p) function of dielectric constant and density of fluid 
G defined by equation (9) 

L length of capillary (m) 

N number of drops/second leaving fluid surface 
P permittivity of fluid (f/m) 

Q/M charge-to-mass ratio of droplet (coul/kg) 

R radius of droplet (m) 

R . radius of perturbation 
V voltage (volts) 
v velocity (m/sec) 

V 0 volume 

AP a pressure rise or drop across meniscus 
€ fluid dielectric constant 
P fluid density (kg/m 3 ) 
cr surface tension (dynes/cm) 


co flow (kg/sec) 
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PARAMETRIC ANALYSIS OF ELECTROSTATIC 
DISPERSION OF A LIQUID 


INTRODUCTION 

The electrostatic dispersion of a fluid is a complicated physical process 
for which no simplified theoretical model has been presented. This analysis 
is an attempt to present a theoretical model which can be used to predict droplet 
radius, charge-to-mass ratio, and velocity, and to identify significant flow- 
dispersion parameters. The results are particularly applicable to analysis and 
evaluation of colloid rocket performance. 

The approach employed in this analysis is to examine the meniscus (liquid 
surface-environment interface) and to devise a model for droplet growth and ex- 
pulsion, with appropriate equations. Then, experimental data are compared with 
the theoretical predictions to validate the model. 


THEORY 

General Description of Flow 

Figure 1 illustrates the overall flow process in the colloid-rocket capillary 
where is the liquid feed pressure, P 2 is the pressure on the liquid side of 
the meniscus, P 3 is the environmental pressure, and V is the applied voltage. 

The pressure drop in the capillary is related to the flow by the Hagen- 
Poiseuille relationship: 


Pj - P 2 = f LA>ci) 2 /2 p A 2 


(1) 


where f = 64/Re for laminar flow. The pressure drop or rise across the meniscus 
is fixed by the La Place-Young equation: 





( 2 ) 


1 



V 



V 


0 - MENISCUS FOR P 3 > P 2 

(2) - MENISCUS FOR P 3 = P 2 • W, > W 2 > W 3 

(3) - MENISCUS FOR P 3 <P 2 

i 

Figure 1. Overall Meniscus Shapes in Flow Capillary 


where R cv is the radius of curvature at the interface. If the meniscus is con- 
vex with c respect to the liquid (case 3 in Figure 1), the liquid pressure P 2 is 
greater than the environmental pressure; if the meniscus is concave with respect 
to the liquid, the liquid pressure is less than the environmental pressure. By 
using equations (1) and (2), the operation of the needle can be visualized. Sup- 
pose Pj is raised to a point where the pressure forces balance the surface 
tension forces, just before the onset of flow: The meniscus would then assume 
position 3 in Figure 1 withPj = P 2 . If a small voltage is applied, resulting in 
a flow, P 2 falls below P x (equation (1)) and closer to P 3 . As the voltage and flow 
are further increased, P 2 will decrease below P 3 , resulting in position 1 in 
Figure 1 (equation (2)). If the voltage is increased even more, the meniscus 
will reach a position from which it can recede no more (because of contact- 
angle and surface-energy effects). After this point, the flow will not increase 


2 


! 


with voltage, and hence is fixed with P x . A colloid-rocket propulsion system 
would ordinarily operate in this regime. 

Analysis 

Figure 2 shows a perturbation, or microscopic fluctuation, on the meniscus 
of the fluid in the capillary, where R £ is the initial radius of the assumed hemi- 
spherical fluctuation. 

The condition for growth of the perturbation is 


Pe > AP ct 


( 3 ) 


where AP^ is defined by equation (2) for R 4 , and P e is the electrostatic pressure 
in the perturbation. 

Assuming the hemispherical perturbation, the electrostatic pressure can 
be expressed by 


Pe 


V 2 

7.4 x 10 5 ttR . 2 


f(e, P) tl] 


( 4 ) 



Figure 2. Perturbation of Meniscus 
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with the following assumptions: 

• The specific electrical conductivity of the fluid is less than 10“ 5 mho /cm. 

• The dielectric constant of the liquid perturbation in the transition surface 
layer changes continuously from a value of E in the interior to 1 in air or 
vacuum, and the transition thickness is about equal to the liquid dipole 
length. 

The factor f( e, p ) can be determined for many fluids by the methods of 
Henriquez [2] and Onsager [3l. 

Substituting equations (4) and (2) into (3) yields 


V 2 

7.4 x 10 s ttR. 2 


f(e. P) 



(5) 


for perturbation growth. Equation (5) may be rewritten as 


R. < 1 

1 14.8 x 10 s TO- 


H e < P ) • 


(6) 


Equation (6) states a condition for growth of a perturbation; i.e., all perturba- 
tions of Rj less than V 2 /(14.8 x 10 s 7 tct f(e, p)) will grow, while all perturbations 
greater than this factor will not. Figure 3 illustrates a distribution curve (type 
unknown) for the perturbations; at voltage V,, all perturbations less than R i maxj 
will grow, and at V 2 all perturbations less than Ri max 2 will grow. For each 
case, the total number of perturbations per second that will grow is expressed 
by 


N 


I 


R.max 


N p dRi 


(7) 


and hence N = area under the distribution curve to the specified limits. 
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V, 1 


V 2 ‘ 



Figure 3. Perturbation Distribution Curve 


From equations (6) and (7), and Figure 3: 


V 2 

NaG — f ( e, p) 


( 8 ) 


where G is a function of the particular perturbation distribution for the meniscus 
of concern. 

Because the distribution curve is unknown, as a first-order approximation 
it will be assumed that G is constant for a given fluid. This in effect assumes a 
flat distribution within the limits of interest. Hence, 


V 2 . 

N = G — f( e, p ) 


(9) 
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Assuming the droplets leaving the liquid surface are of equal size (if they are 
not, then we are considering an average drop), 


co = m 

= N(4/3 77 R 3 p) 


( 10 ) 


for a spherical drop. 

Combining equations (10) and (9) and solving for R yields 


R = 



( 11 ) 


where 


k - G/o- f(e, p) 


The exact mechanism for the expulsion of a droplet from the surface of the 
fluid is unknown. However, Figure 4 is a simplified model of the process which 
consists of a perturbation growing from R A toR, with breakoff atR. At breakoff, 
half the droplet is formed and an equal mass of fluid is extracted from the fluid 
bulk. 

With this model, the charge on the droplet at breakoff can be calculated by 

Q = (2t7R 2 )P(V/R) 

= 2t7 R VP (12) 


where P is the permittivity of the fluid. 
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Figure 4. Model for Drop Breakoff 


The charge-to-mass ratio of a droplet is 


Q_ _ 2RVP 
4/3 v R 3 p 


3 VP 

2 7rR 2 p 


Substitution of equation (11) into (13) yields 

Q U 7/3 

W ~ K i • 2/3 

CO 


(13) 


(14) 


where 


K, = 1.82 


PK 


2/3 


p 1/3 77 1/3 
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The energy equation for the liquid droplet, assuming small viscous dissipa- 
tion and droplet surface energy, is 


QV - "2 M v 2 

Hence, 



(15) 


(16) 


Combining equations (16) and (14) yield 


v 



(17) 


Equations (9) through (17) show that for a known G the significant flow param- 
eters (charge-to-mass ratio, droplet size, velocity, etc.) can be calculated. Thus, 
within the bounds of the theoretical model, G is a parameter of prime significance. 


DATA CORRELATION 

Little data are available to fulfill the initial assumption of a conductivity 
less than l(f 5 mho /cm. The best data found were from twenty-three tests run 
with sodium hydroxide-doped glycerol (resistivity = 7215 ohm-cm) from 
reference (6). 

The parameters Q/M, V, and go were measured during these tests, and G 
could be calculated from equation (14) for each point. Figure 5 shows the re- 
sults of these calculations, which indicate that G is not a function of flow (Pj) 
and shows no discernable variation with voltage. Hence, the assumption of con- 
stant G is validated, although there is some scatter in these data. From these 
data, G is taken as a constant equal to 1.2. 

Using this value for G, the charge-to-mass ratio and velocity of the particles 
were calculated from equations (14) and (17). Figures 6 and 7 show a compari- 
son between these results and the data. 


8 



2.4 


2.2 - 
2.0 - 
1.8 - 


G DEFINED BY EQ. (9) 


O P, = 5" HG 
VP, - lO" HG 
A P | = 20" HG 
ft P , = 30" HG 


DATA FROM REF. [7] 
SODIUM HYDROXIDE - 
DOPED GLYCEROL 


o 



X 1.2 - 
O 

1.0 ~ 


0.8 - 
0.6 
0.4 r- 


0.2 h 


ft 

O 


V 

A 




O 


V 

A 


8 


i 


o 1 —I 1 1 

6 7 8 9 10 


VOLTAGE (kilovolts) 

Figure 5. Variation of Parameter G with Voltage 


Figure 6 validates equation (14): i.e., Q/M varies linearly with V 7/3 /o> 2/3 . 
The assumed model, with constant G, appears to be correct at least within the 
range of these tests (6 to 10 kv, 0.664 x 10” 8 to 4.4 x 10” 8 kg/sec). 

Figure 7 shows a good correlation with the energy equation. The decrement 
in the velocity of these data, as compared to theory, arises from slight errors 
in the reduction of the experimental data. In fact, because the measured charge- 
to-mass ratio is deduced from the energy equation by measuring the velocity, 
the agreement between theory and data in Figure 7 should be as good as in 
Figure 6. 

Reference 6 also contains data from tests run with a sodium ethylate-doped 
glycerol solution. The resistivity of this solution, 3470 ohm-cm, is further 
away from the initial assumption of a greater than 10 s ohm-cm resistivity. 
Nevertheless, Figure 8 shows that equation (14) is valid for this case also. 

If, now, the slope of the data in Figure 8 (Kj) is put into equation (17), the 
velocity of the particles can be calculated. Figure 9 shows the results, in which 
the difference between theory and data again arises from data-reduction errors. 
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1500 


1400 - 
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Figure 6. Charge-to-Mass Ratio vs Flow-Voltage Parameter-Fluid 1 
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Figure 10 is a plot of equation (11) showing the effect of flow and voltage 
on droplet radius for the sodium hydroxide -doped glycerol tests from reference 6. 


DISCUSSION 

Comparison With Other Analysis 


In the present analysis, equation (9) is the criterion used to determine 
droplet radius for a given flow. This method appears more sound than using 
droplet energy minimization (references 5 and 8) which does not consider the 
process or energy expended during formation. For example, Pfeifer, in refer- 
ence 9, computes the charge-to-mass ratio for a droplet by combining an equa- 
tion for the charge buildup on a growing fluctuation with an energy minimization 
criteria on the final droplet. The resulting voltage-flow parameter for the 
calculation of Q/M (to be compared with equation (14)) was (V/o>) 3/7 . Figure 11 
shows data from reference 7 plotted against this parameter. Comparison of 
Figures 11 and 6 shows that the model outlined in this report results in a sig- 
nificant improvement in this correlation 



Figure 10. Droplet Radius vs Flow-Voltage Parameter 
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Figure 11. Charge-to-Mass Ratio vs Flow-Voltage Parameter, Computed According to Pfeifer 
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Effect of Surface Tension 


The parametric grouping of flow and voltage result from the assumed model 
for droplet formation and the constancy of G in equation (9). Identification of 
these parameters provides a realistic basis for comparing test data, as well as 
for reducing the number of required data points. 

The factor G, as stated, is a function of the particular perturbation distribu- 
tion across the surface of the fluid. Thus, G should be dependent upon fluid- 
surface tension. 

\ 

The dependence of G on surface tension can be estimated by assuming that 
the perturbation mechanism is approximated by a small wave across the 
meniscus. The velocity c of the wave is related by 



where R i is the radius of the wave crest, or perturbation. 


(18) 
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The energy of the moving crest is 


E 


o" mv 2 


3V 0 


/R ; 


(19) 


where m is the mass of the crest and V 0 is its volume. 

Substituting expressions for the crest mass and volume in terms of fluid 
density and crest width, w, the energy is expressed by 


E 


1 wR i 


WCT 


( 20 ) 


Hence, for a given energy input, the crest radius is related to surface tension by 


Ri ~ l/cr . 


(21) 


Combining equations (21) and (18) yields 


C ~ cr/ t/p . 


( 22 ) 


The effect of a change in surface tension on N, the number of droplet growth 
sites per second, is: For wave radii meeting the requirement for perturbation 
growth set by equation (6), N should vary directly and in proportion to the wave 
frequency. Hence, 


f 



(23) 


The wavelength, \, is directly proportionaljto R; thus, combining (21) with (23): 


f ~ cr 2 . 


(24) 
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The factor K is proportional to a 2 (with G proportional to a 3 ) and, by the definition 
of Kj in equation (14), 


Q/M ~ a 4/3 . 


(25) 


From equation (21) an increase of surface tension reduces the wave radii. 
Hence, some radii above the critical value (equation (6)) will be reduced enough 
so the perturbations will meet the growth requirements. Hence, Q/M will in- 
crease by a greater amount than indicated by equation (25). 

In summary, the effect of surface tension on Q/M for a given voltage-flow 
parameter can be approximated by equation (25). Again, little data are available 
for correlation. However, the data in Figure 8 show an exponential dependence 
on Q/M on cr of about 1.1. Thus, the analytical estimation of the exponential 
dependence is high when compared with these data. 

Variations in Electric Field and Q/M 


In this analysis, variations in electric -field intensity or droplet charge- 
to-mass ratio across the liquid surface are not considered; average values are 
always considered. 

For the data correlated in Figure 6, the beam efficiency varied from 69.9 
to 95 percent; for the data in Figure 8, the beam-efficiency range was from 
70.3 to 80.2 percent. Within these ranges, the assumption of an average charge- 
to-mass ratio appears to be good. 

Fluid Electrical Conductivity 

An initial assumption in the analysis was that the fluid conductivity is less 
than 10“ 5 mho/cm. Nevertheless, the derived flow-voltage parameter was 
shown to be valid in the 10” 4 mho /cm range. 

\ 

Changing the fluid electrical conductivity does not appear to change the 
flow -voltage parameter, but raises or lowers the floating constant Kj. 

K j also probably depends on the shape of the capillary lip and the fluid- 
metal contact angle. Thus, the difficulty in obtaining a general solution for 
Q/M is apparent. 
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CONCLUSIONS 


The analysis described here has validated the theoretical model, at least 
within the ranges of the data presented. The model serves to identify parametric 
groupings of voltage and flow to determine Q/M, v, and R. These parameters, as well 
as the factor G, provide a basis for comparing colloid dispersion data, for re- 
ducing number of data points, and for predicting performance within groups of tests. 
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